Loss of coherence in double— slit diffraction experiments 
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I. INTRODUCTION 

Coherence is an important building block in physics, 
common to optics and quantum mechanics. In classical 
optics, it simply implies a phase condition related to in- 
terference. Its lack is related to the nature of the source 
of light (non-monocromaticity) and/or the geometrical 
shape of the objects placed on the source-detector path- 
way [H 0- m this sense, incoherence can be regarded 
as a "non dynamical" effect. On the other hand, the 
same concept acquires a more general meaning in quan- 
tum mechanics, where it has to reconcile with the statis- 
tical interpretation of the wave function. When two or 
more states are coherently superposed their properties 
qualitatively differ from those exhibited by the isolated 
components, this being a distinctive feature of the quan- 
tum behavior Q. 

Due to its remarkable relevance in the modern quan- 
tum theory of information and quantum computation Q , 
a topic of much recent interest is the understanding of 
how quantum systems lose their coherence, being deco- 
herence the most widely accepted mechanism [5]. By 
decoherence we understand the irreversible emergence of 
classical properties in a quantum system through its "dy- 
namical" interaction with an environment or bath [(J • 
This idea is better understood by splitting the Hamilto- 
nian describing the full system in three terms: 



H = H A + H B + V AB , 



(1) 



where Ha and Hb are the Hamiltonians governing the 
evolution of system and environment, respectively, and 
Vab accounts for the coupling between them. This last 
term leads to a fast cancellation of the off-diagonal el- 
ements of the system reduced density matrix, which is 
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studied by averaging (i.e., tracing) the full-system den- 
sity matrix over the environment degrees of freedom. 

In general, the visibility of an experimental diffraction 
pattern can be considered the result of combining in- 
coherence and decoherence. In this paper we analyze 
the influence of these two mechanisms using a theoret- 
ical model applied to the double slit experiment with 
cold neutrons carried out by Zeilinger et al. By 
means of a detailed optical study of the neutron beam 
0, we conclude that incoherence and decoherence are 
both needed in order to explain the loss of coherence 
found in the experiment, this being the main result of 
our work. Decoherence is introduced here using a simple 
phenomenological theoretical model that assumes an ex- 
ponential damping of the interferences 9] , and that only 
has one empirical parameter (the coherence degree). It 
should be stressed out that this model is based solely on 
first principles (except for the estimation of the coher- 
ence parameter). Other authors have treated the effects 
of incoherence in similar diffraction experiments with He 
atoms 10] and He clusters 11] from an optical point 
of view, and also the effects of decoherence in fullerene 
diffraction [l^. 

The organization of the paper is as follows. To make 
the paper self-contained, we briefly review in Sec. ITU the 
experiment carried out by Zeilinger et al. |||. Section 
IhI is devoted to the optical description of this experi- 
ment, and in Sec. II VI we present our quantum calcula- 
tions, whose results are compared with the experimental 
data. Finally, the main conclusions derived from this 
work are summarized in Sec. Ivl 



II. 



THE EXPERIMENT 



The double-slit used by Zeilinger et al. [8( consisted of 
a highly absorbing boron wire mounted into the opening 
gap of a single-slit. The dimensions of this arrangement 
(hereafter labelled O) are a\-d-a2 = 21.9-104.1-22.5 /im 
(left slit/boron wire/right slit). The wavelength of the 
neutron beam, AdB = 18.45 A (with a bandwidth AA = 
±1.40 A), was selected by means of an entrance single-slit 
(C) of width w = 20 /mi, located at a distance z = 5 m 
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FIG. 1: Experimental results obtained by Zeilinger et al. p| 
for the double-slit diffraction of cold neutrons (solid circles), 
and theoretical intensity patterns obtained by assuming two 
finite-size, identical slits [see Eq. I|15fl ] (solid line). 



before O. After crossing the double-slit and a pathway 
v = 5 m, neutrons are made to pass through a scanning 
single-slit (D) of width wq = 20 /xm before reaching the 
detector (located behind D). 

The reported results are reproduced in Fig. ^ From 
them, a value for the fringe visibility of V eX p = 0.583 is 
obtained. This magnitude is defined 0,0 as: 



V 



-^max -mi 
Imax ~\~ At 



(2) 



being / max and I m i n the intensities corresponding to the 
central maximum and the first minimum beside it, re- 
spectively. This experimental lack of visibility cannot be 
attributed to simple causes, such as non-convergence of 
the intensity pattern at D or the size difference between 
slits. In the first case, D is far enough from O, and 
the intensity pattern is well-defined (it displays almost 
vanishing minima). In the second case, the difference 
between a\ and is too small (~ 3%), and the inco- 
herence effects due to this asymmetry are irrelevant. In 
the following sections the physical factors leading to this 
high value for the loss of coherence in the experiment are 
examined. 



III. 



THE OPTICAL TREATMENT 



Studies on neutron interferometry can be accurately 
carried out by means of classical optics provided that 
energy is low and spin effects are negligible. Since classi- 
cal optics only accounts for incoherence effects, this ap- 
proach is of great help in our case in order to discern 
whether or not decoherence effects are relevant. This 
analysis will also provide the physical conditions that the 
quantum model for the possible decoherence must fulfill. 
Therefore, we first analyze here the optics of the neutron 



beam to establish the relevance of incoherence effects in 
the diffraction pattern. 

The coherence of light can be understood as the ef- 
fect of stationary stochastic processes i.e., random 
fluctuations in the amplitude and phase of the interfer- 
ing waves of light in time. Although the source of that 
stochasticity is irrelevant for the optical analysis, it turns 
out to be crucial in the quantum theory of decoherence 
5j. The correlation between two processes taking place 
at different points and times is given by the (mutual) co- 
herence function [TEl Il6l H?|. Y(x\,X2\t) 1 which is only 
a function of the time difference, r = t% — due to 
stationarity. From it, the power spectrum [la . Il7| or 
cross-spectral density function 



X\ 



x 2 ;u)= J T(xi 1 X2' ) r)e luJT dr 



(3) 



can be defined. This is the optical analog of the quan- 
tum density matrix. Like it, S(xi,X2',u>) — S* {x2,x\,uj) 
for the off-diagonal elements, while the diagonal ones, 
S(x,x;uj), are real and give the intensity distribution of 
light at a frequency ui. Therefore, the power spectrum is 
the essential element in this optical treatment. 

Since Eq. © refers to the neutron beam, uj can be 
substituted by A without loss of generality. As in Ref . , 
in the theoretical description of the experiment we will 
take into account the pass of the neutron beam through 
C, O, and D. If C is identified with a neutron source 
that can be assumed as incoherent 0, i.e., the intensity 
at each point along the aperture does not depend on any 
other point, then 



S c (xi,x 2 ; A) = q(x 1 ;X)S{xi 



X2 



(4) 



When 



Here, q(x, A) > is the average intensity at ^. wucn 
no absorption takes place at the borders of C (i.e., each 
point of the aperture produces the same intensity) and A 
is independent on the position along this aperture, q(x; A) 
is given by 



q(x;X) 



<A), \x\<w/2, 



(5) 



with s(A) describing the spectral profile of the neutron 
beam. It should be emphasized that to keep the explicit 
dependence on A is important, since the source is not 
totally monochromatic. As will be seen below, the spec- 
tral profile plays a key role on the shape of the measured 
intensity at D within this optical treatment. 

From C to O the neutron beam propagates freely ac- 
cording to the following point-spread function : 



Jk(x-£f/2z 



(6) 



where x and £ refer to the spatial coordinates at O and C 
|19|| . respectively, and k = 2n/X. In this way, the power 
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spectrum just before O becomes 
S' (x 1 ,x 2 ;X) = 



J hco (xi , £1) h co ( x 2 , 6 ) <Sc (6 , 6 ; A)cZ£i d£ 2 

= g(Ax;A) e i *{ x i- a ^V2* j ( 7 ) 



where 



g(Ax; A) = sinc[fc(xi — x 2 )iu/2z]s(A), 



(8) 



with Ax = x\ — X2 and sinc/3 = sin /?//?. As can be seen, 
the finite size of the source gives rise to a coherence in 
the neutron beam, which manifests as a typical single- 
slit diffraction pattern at O. The first minima of this 
correlation pattern |2£j take place at Ax = ±zXdB/w ~ 
±461 /xm. Therefore, the spatial width illuminated by 
its central maximum will be relatively larger than the 
extension covered by the double-slit arrangement. This 
means that the two diffracted beams will display a slight 
divergence, a fact that must be included as a source of 
incoherence in the quantum model (see Sec. II V|) . 

The double-slit acts as a modulator of the incoming 
wave with a modulation function m(x). Thus, the power 
spectrum at O can be written as 

S , o(xi,x 2 ; A) = TO(xi)m*(x 2 )5o(xi,x 2 ; A). (9) 

After crossing the double-slit the two diffracted beams 
described by So propagate freely until reaching D. 
There, the power spectrum, Sd(x%, Xi; A), is computed 
by replacing, in Eq. Q, h C o{x,Q and S c (£i , £ 2 ; A) by 
hoo{x,S,) and So(£i, £2! A), respectively. The diagonal 
elements of Sd give the intensity distribution, I(x;X), 
right before the neutron beam crosses D. For example, 
in the case of two identical very narrow slits (such that 
the modulation function can be expressed as the sum of 
two 5-functions) separated from each other a distance 
d = d + (ax +a2)/2, one obtains 



Sd(%i, Xi] A) = cos(fcAd/2n)s(A) 
+ q(d; A) cos[fc(xi + x 2 )J/2u] 
from which 



^ik(x 1 — x 2 ) /2v 



(10) 



I(x; A) = 1 + sinc(kdw/2z) coa(kdx/v) s(A). (11) 

The separation between two consecutive maxima that re- 
sults from ||TJ is vXdB/d — 73 /im, which is in good 
agreement with the experimental value of Fig. ^ How- 
ever, the value of the fringe visibility, V = 0.881, is still 
far from the experimental value. Moreover, the ampli- 
tude of the oscillations is constant, in strong disagree- 
ment with the experimental pattern. 

If the finite size of D is taken into account, the mea- 
sured intensity for a position x of the detector results 



1 rx+wo/2 

I (x;X) = — / I(x';X)dx', (12) 

w Jx-Wo/2 



which decreases the fringe visibility, but leaves the pro- 
file of the interference pattern unaffected. The non- 
monochromaticity of the neutron beam leads to a sub- 
sequent integration over the spectral bandwidth. Tak- 
ing into account the optical relation k = u>/c, this inte- 
gral can be carried out easily by assuming us = 2ttc/X, 
and substituting the corresponding quantities by AdB and 
AA = 2.80 A. Thus, in the particular case of a uniform 
spectral profile over the bandwidth Aus, one obtains 



Io(x) = / 

A ^ Jus 



Iq(x; u>')du>' 

u-Au/2 

cx 1 + sinc(fcdu;/2z) sinc(fcdu'o/2ii) 
x smc[(AX/X dB )(kdx/2v)] coa(kdx/v), (13) 

where the approximation of constant sinc-functions 
(with A = AdB) has been used, since these functions 
vary slowly within the integration interval (about a 4% 
with respect to their value at A = AdB)- The results ob- 
tained with this expression are much better than those for 
Eq. Hll[) but still far from the experimental data. In par- 
ticular a value of V = 0.772 is obtained, and the averaging 
process introduces a certain modulation with respect to 
the constant-amplitude intensity of Eq. i|ll|) . However, 
the decay in the intensity pattern as the absolute value 
of the detector position increases is not satisfactory. This 
is a general feature [i.e., independent of the spectral pro- 
file determined by s(A)] coming from the assumption of 
point slits at O. 

The modulation of the intensity pattern is a border 
effect caused by the finite size of the slits, while the oscil- 
lations are the result of the interference of two spatially 
separated wavefronts. This can be easily shown by con- 
sidering that the slits are coherently illuminated (i.e., the 
modulation function is described by the sum of two hat 
functions), and their size is relatively small [2l|. In this 
way, one can assume that Eq. (flip still holds, and intro- 
duces the sinc-function caused by a finite-size slit into 
it, thus obtaining 

I{x, A) oc smc 2 (kar]~ /2v) + sinc 2 (fca?7 + /2w) 
+ 2 sinc(fca?y~/2?;) sinc(fca?7~/2i>) 

x sinc(kdw/2z) cos(kdx/v) s(X). (14) 

Here, for the sake of simplicity, both slits are considered 
with the same width, a — (<zi+a 2 )/2, andrj^ = x±vd/2b. 
Using the same arguments, an equivalent expression for 
Eq. 1(13^1 can be calculated, 

Io(x) oc sine 2 (kari~ /2v) + sine 2 (kar] + /2v) 

+ 2 smc(fca77 - /2i>) sinc(fca77 _ /2u) sinc(fcdu>/2z) 

x smc(kdwo/2v) sinc[(AA/AdB)(^d^/2t')] cos(kdxfv), 

(15) 

The corresponding results are shown as solid line in 
Fig. 2] As can be seen, the intensity pattern fits much 
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better the experimental results, presenting only a small 
discrepancy for the central minima. However, the fringe 
visibility, V = 0.760, does not still fully account for the 
value obtained experimentally, although it improves with 
respect to the previous cases. Therefore, the most impor- 
tant conclusion derived from this phcnomenological anal- 
ysis of the experiment is that, although different causes 
contribute to the incoherence of the neutron beam, they 
are not sufficient to explain the experimental results. The 
remaining contribution should then be attributed to dy- 
namical causes, i.e., decoherence, which will be analyzed 
in the next section. 

One final consideration is worth mentioning. The fi- 
nite size of the slits is an important factor to take into 
account in the formulation of the quantum model. In- 
deed, a general phenomenological form for the measured 
intensity accounting for the loss of coherence can be pos- 
tulated from Eq. H15[) as 

I(x) =Ii(x) +l 2 (x) + 2AIi 2 (x)cos(kdx/v), (16) 

where Ti (x) is the intensity function produced by the slit 
i (with i = 1,2), Zi2(x) the intensity function modulating 
the interference, A the coherence degree or damping coef- 
ficient, and the cos-function is the term purely due to the 
interference between the two slits. Notice that coherence 
degree and fringe visibility are not the same unless both 
partial waves contribute equally to the intensity, as hap- 
pens in Eq. (fTTf t . where T\{x) — l 2 (x) = Ii 2 (x) = 1/2. 



IV. THE QUANTUM TREATMENT 

Optical models are attractive because of their their 
simplicity (and accuracy when used in association with 
fitting procedures). However, they do not deal with the 
loss of coherence due to dynamical processes, such as 
those taking place in the regions where the neutron beam 
evolves freely. The effect produced by such processes is 
what we call decoherence. In general, decoherence is un- 
derstood as a complex effect concerning the phase of a 
many-body wave function |22|. For example, according 
to the model proposed by Caldeira and Leggett , de- 
coherence arises from the interaction of the system with a 
bath of harmonic oscillators linearly coupled to the sys- 
tem. The collisions between particles belonging to the 
system and the environment constitute another mech- 
anism proposed to explain the loss of coherence 0, Q 
in physical processes. Although these scattering events 
may be individually quite inefficient, when taking place 
in large numbers they may lead to an exponential damp- 
ing of the off-diagonal elements of the system reduced 
density matrix, since 



dt 



~Ap nm (i), 



(17) 



scattering 



where A = T(l — (^Po\Sl l S n \'ilj )), F is the collision rate, 
and the scattering process between states \n) and m) is 
described by the corresponding S'-matrix element. 



According to the physical considerations made in the 
previous section, in our quantum-mechanical description 
of the process we only consider the evolution of the neu- 
tron beam from O to D, assuming that both outgoing 
partial waves at O are coherent (i.e., finite-size slits). 
Moreover, as mentioned above, these partial waves differ 
in width (because of the different slit sizes), and then 
move apart separating along the two opposite directions 
(notice, however, that the component of this motion is 
relatively small) . These initial conditions are the respon- 
sible for incoherence. Taking this into account, the evo- 
lution of the neutron beam after passing through O (and 
without being affected by an external environment) can 
be described at any subsequent time |24] as 



|tf (0) )t =Ci|Vi)t +c 2 \ip 2 )u 



(18) 



where \ipi}t is the partial wave emerging from slit i, and 
\ c i\ 2 + l c 2| 2 = 1 at any time. In the coordinates rep- 
resentation, the density matrix associated to this wave 
function is 



p| 0) (r,r') = *r(r) ^>(r') 



(o), 



(19) 



with ty\ '(r) = (r|\I> (0 )) t . The diagonal terms of JTHJl give 
the probability density or measured intensity 



(0) (0)/ \ 

Pt = Pt ( r > r ) 



|c 2 | 2 |V'2| 2 + 2|ci||c 2 ||^i |t \*p2 |t cos S t , 

(20) 



where St is the time-dependent phase difference between 
both partial waves. 

When the action of an environment is considered, the 
wave function (|18[) does no longer describe the evolution 
of the full system. Since our main interest focuses on the 
question of how the environment affects the coherence of 
the system by gradually suppressing the (non-classical) 
oscillatory term in Eq. (|20(l . we will assume conditions 
of elastic scattering [25|. Accordingly, only the states 
describing the environment change during the scattering 
process, while the system states remain unchanged. In 
particular, in our case both partial waves will propagate 
as free Gaussian wave packets since they are Gaussians 
initially (see below). In this way, the initial coherent 
state 

|*) = |*(°)> ® \Eo) (21) 

becomes an entangled state with the form 

|*)t = ci|Vi)i ® \Ei)t + c 2 |^2>i ® \E 2 )t (22) 

at any subsequent time. Here, ji/'i) and \ip 2 ) are coherent 
wave packets (or superposition of them), 



x e 



-{z-z ) 2 /4<rl+ip z z 



, (23) 
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where x [z) is the parallel (perpendicular) coordinate 
(with respect to the plane of the double-slit), xq (zo) the 
center of the wave packet along the x (z) axis, p x (p z ) the 
parallel (perpendicular) component of the momentum, 
and a x (a z ) the width of the wave packet along x (z). 
In the case of quasi-plane waves, the corresponding wave 
function can be constructed as a coherent superposition 
of Gaussian wave packets distributed over the width cov- 
ered by each slit. On the other hand, notice from Eq. 12 21) 
that, due to the initial coherence, \E X ) = \E 2 ) = \E ). In 
our case, this initial state for the environment is chosen 
as \Eo) = |I>. 

In order to calculate the diffraction intensity, the re- 
duced density matrix for the system is obtained by trac- 
ing the full-system density matrix over the environment 
states, 



Pi 



J2 t(Ei\pt\Ei)t 



(24) 



In the coordinate representation this operation is equiv- 
alent to carry out the integral of the total density matrix 
over all the 3N degrees of freedom, {ri}fL 1 , of the envi- 
ronment, 



Pt(rJ) = 

J (r,ri,...r„|*) t t (*|r', n, . . . r„) dri • • • dr„. (25) 

By substituting Eq. (|22|l into Eq. I|24(l . one obtains then 
the reduced density matrix: 

Pt(r,r') = 

(1 + M 2 ) [| Cl |Vi t (r)^ t (r') + |c 2 |V 2t (r)^(r')" 
+ 2a tCl c^ lt (r)^ t (r') + 2a t *c 2 c^ 2t (r)V> lt (r'), (26) 

where at — t(-E 2 |i?i)t. From Eq. 126(1 we obtain the 
measured intensity 



Pt = (l + |«t| 2 ) 
|ci| 2 |^i| t 2 + | C2 1 2 1 -02 1 ? + 2A t |ci||c 2 ||V'l|t|V'2|tCOS$ , 



with 



2|a t | 



(1 + M 2 ) 



(27) 



(28) 



For the sake of simplicity, we have assumed that the 
phase difference between the environment states (in- 
cluded in S' t ) is constant. Equation 1(27(1 is the quantum 
analog of the optical (phenomenological) equation 1(16(1 . 
Thus, At is defined as the quantum coherence degree, 
which plays here the same role as A in Eq. i(16() . 

Now, we make the (reasonable) assumption that the 
environment states are too complex to keep mutual co- 
herence as time increases |26| . In this way, even if they 



are initially coherent, they become orthogonal as time 
passes. Thus, one can assume \a t \ — e~ t ^ Tc , being r c 
the coherence time, a measure of how fast the system 
becomes classical (smaller r c implies faster loss of coher- 
ence). By substituting the value of \a t \ in Eq. l(2*5|) one 
obtains 



A t = scch(t/r c ), 



(29) 



which establishes a relationship between coherence de- 
gree and coherence time. Although the value of the co- 
herence time can be derived analytically for interfering 
waves by means of simple Markovian models [9j , Eq. i(29() 
allows to determine it from the coherence degree. In our 
case, we have used A t = 0.63 (or, equivalently, a co- 
herence time of t c ~ 5.08 x 10 -2 s), for which an excel- 
lent agreement between the experimental data and the 
theoretical results is found (see below). This value cor- 
responds to a fringe visibility V = 0.607 (notice that 
A t ~ V, as said at the end of Sec. Illlf) . very close to the 
experimental value (~ 4% above this value). 

Taking all these ingredients into account, we present 
next results for two basic models. In the first one, af- 
ter passing through the corresponding slits, l^il 2 and 
|4' 2 | 2 are considered to be quasi-plane waves; while in 
the second one, their profiles are modelled with Gaus- 
sian functions. We also assume that the slits are on the 
AY"-plane, the direction of propagation of the beam is 
along Z, and £ y £ x (£ x and £ y being the dimensions 
of the slits). In this way, the model can be considered 
as two-dimensional, with the motion taking place in the 
AZ-plane. In order to minimize the effects of the spread- 
ing of the wave function along z with time, given by the 
general expression 



Of = CT 




(30) 



it is enough, in our case, to choose cTq = 2a, what ensures 
a t — (Jq during the time of propagation. 



A. Case I: Quasi— plane waves 

In this first case, the probability density at the left slit 
has been simulated by covering its width with 30 Gaus- 
sian packets, while for the right (wider) one 31 packets 
have been used. For both slits zq = for all Gaus- 
sians, and the distance between neighbor packets has 
been taken equal to a constant, Aj = aj/(ATj — 1), with 
i = 1,2. The width parameter along the A-direction 
is a x ^ = ai/Ni, with Ni = 30 or 31 depending on the 
slit considered, while for the Z-direction cr z = 2a, as ex- 
plained above. The time evolution of the wave function 
has been numerically obtained using Heller's method |27) . 
which is exact in our case. 

In Fig. |2t a) we show the results corresponding to a 
case where only incoherence effects are included [Eq. 1201 ■ 
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section, these oscillations disappear when a model based 
on Gaussians is used. 

To conclude this section, we show in Fig. Hfb) the in- 
tensity resulting when both incoherence and dccohcrcncc 
(with A t = 0.63) effects are included [Eq. (123] ■ As can 
be seen, these results have been greatly improved with 
respect to only consider incoherence, showing an impor- 
tant decrease in the fringe visibility. However, the agree- 
ment with the experimental data is still poor. It should 
be stressed out that they have been obtained by means 
of a reasonable theory based on first principles, with no 
averaging as in Sec. II I II 
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FIG. 2: Intensity obtained from the interference of two quasi- 
plane waves in the case of: (a) only incoherence effects, and 
(b) incoherence plus decoherence effects. For comparison, the 
experimental data (solid circle) are also shown. 



As mention above, after passing through the slits each 
partial wave has a certain motion component along the 
perpendicular direction. The (initial) momentum com- 
ponent along each direction can be estimated using the 
uncertainty principle, according to which p x ~ Ap x ~ 
H/ Ax. Since Ax ~ a, the momentum along the paral- 
lel direction will be p^ — ih/ai, with the plus (minus) 
sign corresponding to the right (left) slit. Hence, the mo- 
mentum along the Z-direction is p£ = [(27rfi,/AdB) 2 — 
(pi^) 2 ] 1 / 2 . As can be seen in the figure, the results in- 
dicate that the effects due to incoherence are very small, 
with the minima of the oscillations almost touching the 
horizontal axis; they display, however, a certain degree of 
asymmetry with respect to x — 0. The problem results 
quite similar to the evolution of the same initial wave 
with both slits of the same width and no perpendicular 
motion. Finally, notice the oscillations that appear for 
x > 400 A. They arise from the diffractive effects caused 
by the borders of the slits, and are not observed in the 
equivalent solid curve of Fig. ^ because of the double av- 
eraging process (with respect to the finite-size of D and 
the bandwidth) carried out there. As will be seen in next 



B. Case II: Gaussian waves 

In this section we show that the results presented in 
the previous section can be highly improved by simply 
considering that the slits are Gaussian [2g, i.e., the wave 
function in Eq. I|18l) is a linear superposition of two Gaus- 
sian wave packets. These Gaussian wave packets are cen- 
tered at Xq 1 ^ = (ai — d)/2 and x^ = (02 + d)/2 (with 
Zq 1 -* = 0), respectively, and their corresponding width 
parameters are taken as <7x = and a z = 2a. The 
value of &X implies that the intensity at the border of 
the corresponding slit is \ipi(±cr x ^ /2)\ 2 /\ipi(xi)\ 2 = c~ 2 

when (x — x{) — u£ /2, i.e., approximately 13.5% of the 
total intensity. Thus, only a very small part of the wave 
will be out of the boundaries defined by the borders of 
the slits. This assumption is in agreement with the ex- 
perimental observations, since Zeilinger et al. reported, 
with regard to the error on the slit widths, that "neu- 
trons penetrating through the boron wire along a chord 
0.2 /xm away from the surface are attenuated by more 
than a factor 1/e", such that these "neutrons would be 
refracted far out of the diffraction pattern" |g| ■ 

In sharp contrast to the case described in the previ- 
ous subsection, when incoherence is considered in this 
model (given, as above, by p^ — ±7i/ a i anc ^ P* = 
l(2irh/\ dB ) 2 -(pi l) ) 2 } 1/2 , for slits with different sizes), the 
resulting intensity, shown in Fig. |3Ja), looks quite differ- 
ent compared with the results of Fig.^a). In this case, it 
can be observed how the Gaussian envelope due to single- 
slit diffraction (which modulates the amplitude of the in- 
tensity maxima) allows the interference pattern to fit the 
shape of the experimental data fairly well. Nonetheless, 
the central minima do not completely agree with the ex- 
perimental ones yet. This theoretical pattern looks very 
much like that given by Eq. (|15l) (see results in Fig. QJ. 
Notice that the diffractive effects have disappeared, and 
then the "tail" of the curve decreases exponentially, as 
expected. Qualitatively, we can state that this model 
reproduces the experimental pattern is as good as the 
optical theory, but with the advantage that here we do 
not need to assume any average. 

Finally, when we consider Eq. I|27|l with A t = 0.63, 
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FIG. 3: Intensity obtained from the interference of two Gaus- 
sian slits for: (a) only incoherence effects, and (b) incoherence 
plus decoherence effects. For comparison, the experimental 
data (solid circle) are also shown. 



or fullerenes in order to explain the behavior of the em- 
pirical data. These models are usually based on reason- 
able physical assumptions about the nature of the par- 
ticle source, and/or the influence of the different objects 
placed in the pathway of the particle beam in the exper- 
imental setup. Here, we have shown how the two-slits 
experiment with cold neutrons performed by Zeilinger et 
al. |8( can be explained to a certain extent by means of a 
simple analytical model based on such kind of consider- 
ations. This analysis has allowed us to characterize the 
relevance of the different sources of incoherence present 
in the experiment. 

However, these procedures mask the presence of an- 
other relevant phenomenon that is likely to exist, and 
may influence importantly the experiment: decoherence. 
This mechanism, also leading to the suppression of coher- 
ence in the system, cannot be described by optical mod- 
els. Decoherence is produced by the dynamical interac- 
tion between system and environment (e.g., by means of 
scattering events), which is not included in such models. 

In this paper we have shown that decoherence is likely 
to exist in Zeilinger et al.'s experiment, and that it can 
explain fairly well with a simple model assuming an 
exponential damping of the interference term Q, and 
using the experimental value of the coherence degree 
(At = 0.63). By means of this model we have been able 
to establish the influence of each element of the experi- 
mental arrangement on the neutron beam. Among these 
influences, the most important one is the effect that the 
double-slit causes on the neutron beam: it splits the ini- 
tial quasi-plane wavefront coming from C into two co- 
herent Gaussian beams. 



which includes both incoherence and decoherence effects 
in this Gaussian-slit approach, the agreement between 
theoretical and experimental results is excellent, as can 
be seen in Fig. EJb). Therefore, this experiment cannot 
only be interpreted in terms of optical models, but needs 
to include decoherence as a mechanism leading to the 
loss of coherence. This is not an a priori expected result 
since the neutron mass is very small. 

V. CONCLUSIONS 

Optical models are currently used in slits diffraction 
experiments with particles like electrons, neutrons, atoms 
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